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ABSTRACT 


We  consider  several  classes  of  interfering  queues  that  appear  in  packet-radio 
networks.  We  analyse  the  class  of  systems  where  one  of  the  queues  is  given  full 
priority  and  obtain  an  expression  for  the  joint  probability  distribution  of  the  queue 
lengths.  For  ALOHA-tvpe  systems  with  two  symmetric  queues  we  calculate  the  average 
packet  waiting  time  and  queue  lengths  and  for  symmetric  systems  with  an  arbitrary 
number  of  subscribers  we  develop  a  method  to  approximate  these  quantities.  The 
approximation  turns  out  to  be  close  to  the  analysis  and  simulation  results. 
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1.  Introduction 


The  present  study  was  motivated  by  the  problem  of  investigating  the  behavior 
of  random-multiple-access  systems  and  of  packet-radio  networks.  These  systems  are 
characterized  by  the  fact  that  a  number  of  radio  stations  exchange  digital  informa¬ 
tion  by  using  a  distributed  random  access  algorithm  on  a  common  radio  channel.  In 
such  situations,  whenever  a  given  station  attempts  transmission  of  a  packet  to 
another  station,  the  attempt  may  be  unsuccessful,  in  which  case  the  packet  must  be 
retransmitted.  In  addition  to  channel  noise,  unsuccessful  transmissions  occur 
because  of  interference  from  another  station  trying  to  send  a  packet  over  the 
common  channel  at  the  same  time  or  by  the  fact  that  the  intended  receiver  is  itself 
in  transmission  mode,  in  which  case  it  is  not  able  to  detect  incoming  packets.  The 
fact  that  the  activity  at  one  node  affects  the  behavior  of  the  queue  at  other  nearby 
stations  gives  rise  to  statistical  dependence  between  the  queues  at  the  nodes. 

Unfortunately,  in  general, the  queue  length  statistical  dependence  is  quite 
complicated  and  there  is  little  hope  to  obtain  explicit  analytical  results  for  general 
topology  networks.  The  purpose  of  this  paper  is  to  present  several  analytic  as  well 
as  approximation  results  for  certain  classes  of  interfering  queues.  We  assume 
throughout  the  paper  that  all  packets  have  equal  length  and  that  the  time  is  divided 
in  slots  corresponding  to  the  transmission  time  of  a  packet.  A  station  may  start 
packet  transmissions  only  at  the  beginning  of  a  slot  and  the  distances  between  stations 
are  assumed  to  be  such  that  propagation  delay  is  negligible.  Also,  we  neglect  channel 
noise  and  assume  no  channel  errors. 

Because  of  the  difficulty  in  the  analysis  of  dependent  queues  of  the  type 
introduced  above,  even  the  case  of  two  queues  cannot  be  treated  analytically  in  the 
general  case.  However,  we  consider  here  two  classes  of  systems  with  two  dependent 
queues  where  such  results  can  be  obtained.  The  first  is  the  case  when  the  length  of 
one  of  the  two  queues  is  not  allowed  to  decrease,  unless  the  other  queue  is  empty.  A 
two  node  ALOHA  system  where  one  of  the  nodes  is  given  full  priority  is  an  example  for 


such  a  situation;  other  examples  are  given  in  Sec.  2.  No  other  restrictions  are 
necessary  in  order  to  allow  for  analytical  solution  of  this  class;  in  particular, 
the  inputs  may  have  arbitrary  distributions  and  need  not  be  independent  processes. 

For  this  class  of  problems  we  present  a  general  methoa  for  deriving  the  generating 
functions  of  the  queue  lengths  and  of  the  average  delay  times.  Then  these  general 
results  are  applied  to  three  special  cases  of  packet-radio  networks  that  can  be  shown 
to  belong  to  the  considered  class  of  systems. 

The  second  class  of  problems  for  which  we  can  obtain  explicit  analytical 
results  is  the  case  of  a  two  node  symmetric  ALOHA  system.  For  this  situation  we 
cannot  obtain  the  queue  length  probability  distribution  (or  generating  function), 
but  we  give  a  method  for  calculating  the  average  queue  length  and  hence  the  average 
time  delay.  The  results  are  given  in  Sec.  3. 

Since,  as  said  before,  explicit  analytical  results  are  hard  to  obtain  for  more 
general  situations,  another  way  to  approach  the  problem  is  to  obtain  good  approxima¬ 
tions.  An  approximation  method  applicable  to  a  symmetric  ALOHA-type  system  with 
arbitrary  number  of  stations  is  introduced  in  Sec.  4.  We  obtain  there  the  approximat 
average  queueing  delay  in  such  systems  and  compare  this  with  the  exact  result  for  two 
node  networks  and  with  simulation  results  for  larger  networks. 

Discrete  time  systems  involving  interfering  queues  have  been  rarely  treated  in 
the  literature.  In  [2]  a  loop  system  using  Asynchronous  Time-Division  Multiplexing 
has  been  analysed.  In  this  system  user  i  may  use  a  slot  for  transmission  only  if 
all  users  l,2,...,i-l  have  nothing  to  transmit.  In  [3],  the  author  investigates 
the  case  of  two  queues  in  tandem  where  each  queue  always  attempts  transmission 
provided  it  has  a  message  in  the  buffer.  Both  systems  considered  in  [2]  (for  two 
users)  and  [3]  are  cases  that  belong  to  the  first  general  class  of  systems  considered 
in  the  present  paper.  In  [4]  another  system  of  two  interfering  queues  is  considered, 
whereby  only  one  event  (i.e.  an  arrival  or  departure)  may  occur  during  a  given  slot. 


1 

I 


This  system  has  been  shown  to  have  the  product  form  solution.  Finally,  we  may 
mention  [5]  and  [6],  where  a  slotted  ALOHA  network  with  finite  number  of  users  has 
been  examined  and  a  method  was  suggested  for  obtaining  an  approximate  solution  for 


mis  svstem. 
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2 .  Analytical  Results 

In  this  section  we  consider  a  class  of  discrete-time  queueing  systems  consisting 
of  two  queues  with  the  following  properties  :  packets  arrive  randomly  at  the  queues 
from  two  sources,  that  in  general  may  be  correlated.  Let  A^(t)  and  A2 Ct )  be  the 
number  of  packets  entering  node  1  and  node  2  from  their  corresponding  sources  in 
the  time  interval  (t,  t+1] .  The  input  process  [(A  (t),  A2 Ct ) ]  is  assumed  to  be  a 
sequence  of  independent  identically  distributed  random  vectors  with  integer-valued 
elements.  Let 

CO  00 

a(i, j)  =  Prob(A  (t)=i ,A?(t)=j)  ;  £  £  a(i,j)  =  1  .  (1) 

1  i=0  j=0 

and 

F(x,y)  =  E[xAl(t)yA2(t)]  »  I  [  a(i,j)xiyi  (2) 

i=0  j«0 

We  assume  that  F(x,v)  cannot  be  x - independent ,  namely  that  packets  arrive  at  the 
first  queue  with  nonzero  probability  and  that  the  queues  have  infinite  buffers. 

Next,  we  describe  the  departure  processes.  It  is  assumed  that  no  more  than 
one  packet  may  leave  each  queue  in  any  given  time  slot  and  the  combined  departure 
process  is  taken  to  be  as  follows  :  When  both  queues  are  empty,  no  departures  may 
occur  (packets  arriving  during  a  given  slot  may  depart  only  in  the  next  one) .  When 
only  one  of  the  queues  is  nonempty,  a  departure  from  that  queue  may  occur  and  the 
packet  may  be  transferred  either  to  the  outside  of  the  system  or  to  the  other  queue. 

We  denote  by  p^,  p^  the  respective  probabilities  when  the  nonempty  queue  is  queue 
2  and  by  P^q.  pJq  the  corresponding  probabilities  when  the  nonempty  queue  is  queue  1. 
Clearly  1  -  p^  -  p^  is  the  probability  that  no  departure  occurs  from  the  nonempty 
queue  2.  The  specific  class  of  dependent  queues  considered  in  this  section  is 
characterized  by  an  assumption  on  departures  when  both  queues  are  nonempty.  For  this 
case  it  is  assumed  that  a  departure  may  occur  only  from  queue  2.  We  denote  by  p^ , 
pjj  the  probabilities  that  the  departing  packet  leaves  the  system  or  joins  the  other 
queue  (queue  1)  respectively. 
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Consider  the  steady-state  joint  generating  function  of  the  queue  lengths  : 


r L.ct)  L2(tf 

C(x,y)  =  lim  Ex  y  (3) 

t-*-oo  I—  -  ■ 

where  L^Ct),  L 2(t)  are  the  queue  lengths  at  time  t  at  nodes  1  and  2  respectively, 
and  where  we  assume  that  the  Markov  chain  [L^(t),  L^Ct)]  is  ergodic,  namely 
G(i0,0)>0.  For  the  class  under  consideration  we  can  compute  the  function  G(x,y)  and 
in  the  Appendix  it  is  shown  that  G(x,y)  has  the  following  form  : 


rev  v)  -  FTy  vi  •  btXyjGCx.O)  +  c(x,y)G(0,y)  +  d(x,y)G(0,0) 
,>J  ~  1  ,>j  x •  e ( x , y] 


where  the  functions  b(x,y),  c(x,y),  d(x,y) ,  e(x,y),  G(y,0),  G(0,y)  and  the  constant 
G(0,0)  are  defined  in  the  Appendix. 

This  general  form  can  be  ma<Te  more  explicit  for  certain  special  cases.  We 
next  consider  several  examples  of  two  node  packet  networks,  where  it  turns  out  that 
the  general  assumptions,  given  earlier,  characterizing  the  class  of  two  dependent 
queues  indeed  hold.  The  networks  under  consideration  are  given  in  Fig.  1.  In  all 
cases  the  nodes  share  a  common  radio  channel  and  are  equipped  with  radio  transmitter 
devices  and  in  systems  2  and  3,  node  2  has  also  a  receiving  device.  Node  2  can  either 
transmit  or  receive,  but  not  simultaneously.  The  circle  in  Fig.  1  represents  a  station 
equipped  with  a  radio  receiver  that  receives  packets  correctly  provided  that  there  is 
no  interference.  Finally,  instantaneous  feedback  to  the  transmitter  is  assumedi 
meaning  that  a  transmitter  knows  at  the  end  of  the  slot  if  the  packet  has  been  received 
correctly. 

In  all  three  systems  of  Fig.  1  node  2  is  assumed  to  have  full  access  capability 
to  the  common  channel.  This  means  that  it  always  transmit  a  packet  when  its  buffer  is 
not  empty,  while  if  its  buffer  is  empty  the  node  does  not  transmit  and  in  systems  2 
and  3  it  is  able  to  receive  packets  transmitted  by  node  1.  Node  1  has  only  partial 
access  capability  to  the  channel  and  its  transmission  policy  is  randomized  as  follows  : 
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At  the  beginning  of  each  slot  for  which  its  own  buffer  is  nonempty,  node  1 
tosses  a  coin  with  probability  of  success  p,  independently  of  any  other  event 
in  the  system  and  in  case  of  success  the  node  attempts  to  transmit  the  packet  at 
the  head  of  the  queue.  Both  nodes  are  able  to  detect  at  the  end  of  the  slot  if 
their  transmissions  were  successful.  At  any  node,  if  the  transmission  is  not 
successful,  either  because  the  packet  was  sent  to  the  other  node  while  that  node  was 
not  ready  to  receive  it,  or  because  of  interference  with  a  packet  transmitted  by  the 
other  node,  the  transmitter  repeats  the  procedure  described  above. 

Since  node  2  has  full  access  capability  to  the  channel  and  node  2  cannot 
receive  and  transmit  packets  at  the  same  time,  it  is  clear  that  in  all  cases  no 
packets  can  leave  node  1  whenever  the  queue  at  node  2  is  nonempty  and  therefore  all 


cases  of  Fig.  1  belong  to  the  class  of  queues  considered  earlier  in  this  section. 

k 

We  now  turn  to  calculate  the  parameters  (p„,  Osi.j.ksl,  i+j>0}  in  each  of 
these  three  systems.  System  1  depicted  in  Fig.  1(a)  represents  a  two  node  non- 


symmetric  ALOHA  network,  where  both  nodes  send  their  packets  to  the  station.  Since 
no  packets  are  sent  from  one  node  to  the  other  we  have  =  pj^  =  0.  When 

one  of  the  nodes  has  packets  to  transmit  while  the  other  is  empty,  any  attempted 


transmission  is  successful.  Since  node  2  transmits  with  probability  1  and  node  1 
with  probability  p,  we  have  p^Q  =  p;  p^  =  1.  When  both  nodes  have  nonempty 
queues,  successful  transmission  occurs  at  node  2  whenever  node  1  does  not  attempt 
transmission  and  therefore1  p^  =  p.  System  2,  depicted  in  Fig.  1(b)  represents  a 
situation  of  two  tandem  nodes  where  the  station  that  is  the  "sink"  for  the  packets 


transmitted  by  node  2,  is  out  of  the  transmission  range  of  node  1.  Therefore  node  1 


cannot  interfere  with  the  transmissions  of  node  2.  However  node  2  does  interfere 


with  the  transmissions  of  node  1  since  when  it  is  transmitting,  it  does  not  accept 
packets  transmitted  by  node  1.  Consequently  pj^  =  p;  p^  =  p^  =  1; 

=  ^01  =  Pli  =  System  3,  depicted  in  Fig.  1(c)  differs  from  system  2  only  in 

that  the  station  is  in  the  transmission  range  of  node  1,  therefore  node  1  does 


7 


interfere  with  the  transmissions  of  node  2  in  this  case,  and  therefore  the  parameters 


ar : 


i  0.0-0  1 

P10  “  p’  P01  =  1>  P11  "  p’  pI0  =  p01 


'll 


=  0. 


Numerical  Results 

Although  the  results  of  this  section  hold  for  general  input  processes,  the 
equations  become  much  simpler  when  one  considers  independent  Bernoulli  processes. 
In  this  case  we  have 


F(x,y)  =  (xrx  +  r^  (yr2  +  r2)  (5) 

where  r ^ ,  r0  are  the  input  rates.  For  this  example  we  calculate  for  each  of  the 
three  systems  the  average  delays  (in  units  of  slots)  T  ,  T2  at  nodes  1  and  2 
respectively  and  the  total  average  delay  T  in  the  network.  This  is  done  by  first 
calculating  the  average  queue  lengths  at  the  nodes  and  then  applying  Little's  Theorem 
[7],  After  straightforward  but  tedious  algebra  the  following  results  are  obtained  : 
System  1  (Fig.  1(a)) 

(P)2+r2P  ri r2PP 

Tj  =  1  +  - - — r  +  - j-1 -  C6) 

p(p-r2)-r1p  (p-r2)  [p(p-r2) -r^] 


T2  =  !  + 


‘lp 


and 


(P-r2)“ 
rlTl  +  r2T2 


rl+r2 


(7) 


(8) 


where  these  equations  hold  for  p(p-r2)  >  r^p  which  is  the  ergodicity  condition  in 
this  system. 

In  Fig.  2,3,4  we  plot  T^,T2  and  T  respectively  versus  p,  the  transmission 
probability  at  node  1,  for  r^  =0.1  and  r2  ranges  from  0.01  to  0.4.  As  expected 
the  average  delay  at  node  2  increases  when  p  increases  since  then  its  transmissions 
are  more  frequently  interfered  with  transmissions  from  node  1.  More  interesting  is 
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the  behavior  of  the  average  delay  at  node  1.  Here  p  has  some  value  for  which 

is  minimized  (for  given  r^  and  r7)4  When  p  either  increases  or  decreases  from  this 

value  increases.  The  reason  is  that  when  p  becomes  small,  node  1  attempts 

to  transmit  relatively  rarely,  so  its  queue  increases.  When  p  becomes  large,  then 

node  1  attempts  to  transmit  more  frequently,  thus  interfering  with  the  transmissions 

of  node  2,  and  the  queue  lengths  at  both  nodes  are  large.  As  we  see  from  Fig.  4  the 

parameter  p  is  a  very  critical  design  parameter  of  this  system  and  for  given  values 

for  r^  and  r,,  there  exists  an  optimal  p,  denoted  by  p*,  that  minimizes  the 

total  average  delay  in  the  network.  In  Fig.  5  p*  is  plotted  versus  r^  for  various 

values  of  r7.  Notice  that  p*  is  much  less  sensitive  to  changes  in  r^  than  to 

changes  in  r7.  When  r7  is  small,  then  p*  should  be  large  for  all  values  of  r  , 

since  interference  between  transmissions  is  rare.  When  r2  increases,  p*  should 

decrease  in  order  to  reduce  the  interference.  Finally,  in  Fig.  6,  T  .  the  minimum 

mm  - 

total  average  delay  is  plotted  versus  y  the  total  throughput  of  the  system,  when 
r^  =  r7  =  r  (clearly  y  =  2r) . 

System  2  (Fig.  1(b)) 

The  average  delays  are  : 


T1  =  1  + 


rlP+f2(1"p?2^ 

r2[p(l-r1-r2)-ri; 


T7  =  — - — 
2  r1+r2 


rT  +  I - 

2  1-r¬ 


and 


T  =  — i —  T.  +  T, 
rl+r2  1  2 


(9) 


(10) 


(ID 


where  these  equations  hold  for  p(l-r1-r2)  >  r^  which  is  the  ergodicity  condition 
for  this  system. 

In  this  system  node  1  does  not  interfere  with  the  transmissions  of  node  2. 
Therefore,  it  is  optimal  to  always  attempt  transmission  at  node  1  as  well,  namely  to  , 

,Cf 

take  p  -  1.  Eq.  (9)  indeed  shows  that  T^  is  monotonically  decreasing  when  p  y 
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increases,  and  it  achieves  its  minimum  value  for  p  ■  1.  In  Fig.  7,8  and  T  are 

plotted  respectively  versus  r^  for  p  =  1  when  r2  ranges  from  0.01  to  0.8.  T2 
is  not  plotted  here  since  it  takes  values  from  1  to  1.5  only  when  0.01^r2s0.8, 

0  <  r,  <  0.5.  From  Fig.  7  we  see,  as  expected,  that  Tj  increases  when  the  arrival 
rate  r^  to  node  1  increases.  When  r0  increases,  then  the  time  spent  by  node  2 
in  the  transmitting  mode  increases.  Therefore,  packets  transmitted  by  node  1  are 
rarely  received  by  node  2  and  the  average  time  delay  T  at  node  1,  increases  in 
this  case  as  seen  in  Fig.  7. 

System  5  (Fig.  1(c)) 

pir.+r +  (p-r?)(p)2  pt^+r^r^r^)  ] 

T .  =  1  *  - - — — - - - - 5 - I12) 

[p(p-r,)-ri](p-r2)  [?1f2)“-p(l-r1r2)](p;r2) 

f 

1  |Vr2f2  pCr1^r2(r1+r1r2)][p(p-r2)-r1]  ^ 

2  rl+r-  p-r2  (p-r2)  [r1(r2)2-p(l  -  r^)] 


and 
T  = 


r, 

i. 

- T, 

rl+r2  1 


To 


(14) 


where  these  equations  hold  for  p(p-r2)  -  r^  >  0  which  is  the  ergodicity  condition 
for  this  system.  Similarily  to  System  1,  T ^ ,  T2  and  T  are  plotted  in  Fig.  9,10,11 
for  System  5.  The  behavior  of  the  average  delay  in  this  system  is  very  similar  to 
that  of  System  1,  therefore  we  will  not  give  here  more  detailed  explanations. 
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5 .  Symmetric  Two-Node  ALOHA  Network  Delay  Analysis 

In  this  section  a  two  node  symmetric  ALOHA  network  is  considered.  This  network 
is  similar  to  System  1  with  the  following  modifications.  Here  node  2  uses  the  same 
channel  access  scheme  as  node  1,  i.e.  at  the  beginning  of  each  slot,  if  its  buffer  is 
not  emptVj  node  2  tosses  a  coin,  independently  from  any  other  event  in  the  system, 
with  probability  of  success  p.  According  to  the  outcome,  the  node  either  transmits 
or  remains  silent  during  the  current  slot.  In  addition,  it  is  assumed  that  the 
arrival  processes  to  the  nodes  are  independent  Bernoulli  processes  with  equal  rates, 
denoted  by0<r<l/4.  Therefore  we  have  for  this  system  F(x,y)  =  (xr+r) (vr+r) .  By 
the  same  method  used  in  the  Appendix  it  is  easy  to  see  that  for  the  current  system 
we  have  : 


G(x,y)=F(x  y)p[xl-'y-1)-p(2xy-x-y)lG(x,.Q)  +  [yCx-l)-p(2xy-x-y)]G(0,y)^pf2xy-x-y]G(0,0) 

xv-F(x,y) [ (x+y)pp+xy (p'+ (p) ") ] 

In  this  case  we  cannot  obtain  an  explicit  form  for  G(x,0),  G(0,y),  G(0,0)  and  hence 
for  G(x,y).  However  we  can  exploit  the  symmetry  to  obtain  an  expression  for  the 
average  delay  in  the  system.  If  we  denote  by  Gj(x,y),  G9(x,y)  the  derivative  of 
G(x,y)  with  respect  to  x  and  y  respectively,  we  clearly  have  G^(l,l)  =  G-,(l,l) 
and  G^(1,0)  =  Go(0,l).  Then  from  (15)  we  obtain 


_  -> 


Gjd.l)  =  r 


r[p+(p)~]  -  P  G^ ( 1 ,0) 


and 


PP  -  r 


^[G(x,x)] 


2r  + 


Gj  (1 ,0)p(l-2p)  r“  +  2r-4rpp 


x=l 


pp  -  r 


2 (pp  -  r) 


(16) 


(17) 


for  pp  >  r. 

Now,  if  we  use  the  fact  that 


^[G(x,x)] 


x=l 


=  G1(l,l)  +  G:(l,l)  =  2G1(1,1) 


(18) 


(15) 


we  can  solve  for  G^(1,0)  and  hence  for  G^(l,l)  and  obtain 

,-,2  1 

(P)  +  7  rp 

GjCl.l)  =  G2(l,l)  =  r  +  r - - - - -  (19) 

PP  *  r 

Therefore,  applying  Little's  Theorem  we  obtain  the  average  system  delay 
G,  (1,1)  (p)2  +  j  rp 

T  =  -  =  1  +  -  for  pp  >  r  .  (20) 

pp  -  r 

2  1/2 

From  (20)  it  is  found  that  p*  =  l-(0.5r  +  [0. 5r(l-r+0. Sr  )]  }/(l-0.5r)  minimizes 

T  for  0<r<0.25.  In  Fig.  6  the  minimum  total  average  delay  is  plotted  versus 

Y,  the  total  throughput  of  this  system.  Comparing  the  curves  in  Fig.  6  it  is  clear 
that  the  non-symmetric  access  scheme  used  in  system  1  provides  very  slight  better 
performance  than  the  symmetric  access  scheme,  when  the  arrival  rates  into  the  nodes 
are  equal  (the  difference  in  the  minimum  delay  is  less  than  3.5%  in  the  range  0<y<0.5). 
However  the  former  scheme  is  unfair,  giving  priority  to  queue  2,  although  the  arrival 
rates  into  the  two  queues  are  the  same. 

Finally  we  may  mention  that  the  method  presented  in  this  section  for  calculating 
average  delay  times  without  obtaining  explicitly  the  generating  functions,  can  be  used 
in  many  other  symmetric  (i.e.  G(x,y)  =  G(y,x))  two  node  systems.  Specifically,  we  can 
easily  obtain  average  delay  time  in  the  symmetric  two  nodes  ALOHA  network  for  general 
arrival  processes  into  the  nodes  for  which  F(x,y)  =  F(y,x). 
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4.  Symmetric  ALOHA  Network  with  M  Nodes  :  Approximate  delay 

Consider  an  ALOHA  network  consisting  of  M  nodes,  that  share  a  common  radio 
channel.  Assume  that  all  nodes  in  the  network  use  the  same  channel  access  scheme 
as  the  two  nodes  in  the  previous  section,  i.e.  each  node  tosses  independently  a 
coin  with  the  same  probability  of  success  p  at  the  beginning  of  each  slot  when 
its  queue  is  nonempty.  Also  assume  that  the  arrival  process  at  each  node  is  a 
Bernoulli  process,  independent  from  node  to  node,  with  rate  r.  For  this  case 
we  cannot  obtain  the  exact  average  delay  for  M>2,  and  therefore  we  must  consider 
approximations.  The  approximation  method  proposed  here  is  the  following.  Each  node 
i  in  the  network  may  at  any  time  be  in  one  of  M  possible  situations.  Situation 
j  refers  to  the  case  where  j-1  nodes  other  than  i  have  nonempty  queues,  while 
the  other  M-j  nodes  (not  including  i)  are  empty.  The  approximation  considered  here 
consists  of  assuming  that  in  steady  state,  while  in  situation  j,  node  i  behaves 
as  a  discrete  M/M/1  queue  with  arrival  parameter  r  and  departure  parameter 
p(l-p)J_1.  If  the  transitions  between  the  various  situations  are  neglected  then  the 
average  number  of  packets  at  node  i  denoted  by  L^  is  : 


L.  =  >9.  -  v  .  — 

1  j  =  l  J  p C 1 -p) ^ ' 1 -r 

for  i=l,2,...,M  and  r<p(l-p)M~\  where 
situation  j.  We  approximate  9^  as 


0  . 
J 


(21) 


is  the  probability  of  being  in 


9  . 
J 


M-j  r  j-1 
P 


j  =  1,2,.  ...M 


(22) 


where  —  approximates  the  probability  that  a  node  has  packets  ready  for  transmission 
and  we  assume  independence  between  the  nodes.  Using  (21)  and  (22)  and  applying 


Little's  Theorem  we  find  that  the  approximate  total  average  delay  T  is  given  bv  : 

aP 


T 

ap 


M 

I 

j  =  l 


P 


1  (1  -  I)»H 

P 


1-r 

p(l-p)^  _1-r 


(23) 
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The  general  formula  (23)  can  be  specialized  to  the  case  M=2,  and  for  this  case  we 
can  compare  the  approximate  and  exact  results.  For  M=2  we  have  from  (20) 


and  from  (23) 


T 


analysis 


, -.2  1 
(P)  +  2  ^ 

PP  -  r 


T 

ap 


+ 


(P)2  +  rp 
PP  -  r 


(24) 


(25) 


T  and  T  ,  are  plotted  versus  p  in  Fig.  12,  for  various  values  of  r. 

ap  analysis  r 

It  can  be  seen  that  although  we  have  used  a  simple  approximation,  it  is  quite  close 
to  the  exact  values.  We  also  compare  this  approximation  versus  simulation  results, 
for  networks  having  three  and  four  nodes.  This  comparison  is  plotted  in  Fig,  13 
and  14  and  shows  again  good  results. 
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Footnote; 

1.  We  use  the  notation  v  =  1-v. 


The  subscripts  1  and  2  correspond  to  the  derivative  of  the  function  with 
respect  to  the  first  and  second  variables  respectively. 


In  order  to  describe  the  behavior  of  the  queue  lengths,  we  need  several 
definitions.  Let  D^^t)  be  a  binary-valued  random  variable  that  takes  value 
I  if  L  (t)  =  0,  L2 (t)  >  0  (L^(t)  and  L2 (t )  are  the  queue  lengths  at  time 

t  at  nodes  1  and  2  respectively) ,  a  departure  occurs  (from  queue  2)  and 
the  packet  leaves  the  system.  Similarly  D^(t)  takes  value  1  under  the  same 
conditions,  except  that  the  packet  is  transferred  to  queue  1.  In  a  similar  way 
we  define  the  binary-valued  variables  D^ft),  D^Ct)  for  the  case  that 
L1(t)>0,  L2(t)=0  and  D^(t),  D^(t)  for  the  situation  when  both  queues  are 

nonempty  and  departure  occurs  from  queue  2.  Then  the  recursive  equations  for 

L. (t)  are  : 

1 

Ljft+1)  =  • 


L:(t+1)  =  ■ 

From  (Al),  (A2)  and  (3)  we  have  for  t-*»  that 

G(x,y)  =  F(X,y){G(0,0)+[G(x,0)-G(0,0)][x-1pJ0+x*1yp}0+  (p°Q+p}0)]  + 

+  [G(0,y)  -G(0,0)  ]  [y” ^Q^y" lxPoi  +  (p01+P01)J  + 

+  [G(x,y)-G(x,0)-G(0,y)+G(0,0)]  [y^p^+y'^p^  +  (p^+p^)]}  (A3) 


Aj (t)  if 

A1(t)+L1(t)-D°0(t)-D}0(t)  if 
A1(t)+Dj1(t)  if 

A1(t)+L1(t)+Dj1(t)  if 

A2(t)  if 

A2(t)+D}0(t)  if 

A2(t)+L2(t)-D°1(t)-Dj1(t)  if 

A2(t)+L2(t)-D°1(t)-Dj1(t)  if 


LjCO-O,  L2(t)-0 
L  (t)>0,  L~ (t) =0 

(Al) 

L^ (t)=0,  L2(t)>0 
L1(t)>0,  L2(t)>0 

Lj (t) =0,  L2(t)«0 
(t) >0,  L2(t)=0 

Lj^ (t)  =0,  L2(t)>0  (A2) 

I-,  (t) >0,  L- (t)>0 


(Remember  that  v  denotes  1-v). 
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Arranging  (A5)  we  obtain 


Gfx  v)  =  F(x  y)  b(x>y)G(x»Q)+c(x»y)G(0,y)+d(x,y)G(0,0) 
i.  >yj  *-  x  e(x,y) 


(A4) 


where 


b(x>y) 

c(x,y) 

d(x,y) 

e(x,y) 


=  y(pio+yp!o}  -  x(P?i+xPn)+xy(P?i+Pn-P?0-Pio) 

(AS) 

=  xtPorPn+x(pJrPn)]  +  xy(p?i+Pn-PorpJi) 

(A6) 

,  o  1  ,  r  0  0  .1  1,.  ,0  1  0 

=-y(Pio+yPio)^x[Pii-Poi*x(Pii-p0i)]  +  xy(p10+p10+p01 

10  1, 
+P0l'Pll'Pll^ 

=  y-F(x,y)[pJ1+xpJ1+y(pJ1+pj1)] 

(A8) 

From  (A4)  we  see  that  the  steady-state  generating  functions  for  the  queues'  lengths 
are  : 

For  queue  1  : 


G (x, 1) 


Ffx  n  b (x, 1) G(x,0) +c (x , 1)G(0, 1) +d(x, 1)G (0,0) 
1  J  x  e(x,l) 


For  queue  2  : 


G(l,y) 


FM  vl  b(l,y)G(l,Q)+c(l,y)GC0,y)+dCl,y)G(0,0) 
U,XJ  e(l  ,y) 


(A9) 


(A10) 


We  still  have  to  determine  the  functions  G(x,0),  G(0,y)  and  the  constant  G(0,0). 
To  find  G(0,y)  let  x-K)  in  (A3).  Then 

G(0,y)  =  F(0,y){G(0,0)+G1(0,0)(p°0+ypJ0)  + 

+  [G(0,y)-G(0,0)][y-1p°1+  (Po^P^)]}  (All) 

where  G^(x,y)  is  the  derivative  of  G(x,y)  with  respect  to  x. 


From  (All)  we  obtain 


G(0,y) 


F(0,y) 


1  -  F(0,y) [y_1p§1  +  (pg1+pj1)  ] 


(A12) 


Before  proceeding  we  need  to  prove  the  following  Lemma. 
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Lemma 

For  | x | < 1  the  equation 

e(x,y)  =  0 

where  e(x,y)  is  the  expression  of  (A8)  has  a  unique  solution  within  the  unit 
circle  |y|=l. 

Proof 

Since  F(x,y)  is  not  independent  of  x,  then  there  exists  a(i,j)>0  for  some 
i  and  some  j>0.  Therefore  for  |x|<l  and  jyj  =  1  we  have 

|F(x,y)[p°1+p}1+y  (P^-Pji)]  I  * 

00  00  00  00  OO  CO 

s  i  F(x,y) |  =  11  l  a(i , j)x1y^  j  $  l  l  aCi, j) I x J 1  <  l  l  (a(i, j)  =  l=|y|  (A13) 

i=0  j=0  i=0  j=0  i=0  j=0 

Hence,  applying  Rouches'  theorem  [8],  e(x,y)  has  exactly  one  zero  within  |y|=l  for 
i  x .  <  1. 

Q.E.D. 

Now  let  t  be  the  solution  of 

FCO.tHt*1?^  H-  (pj^pjj)]  =  1  (A14) 

in  the  unit  circle.  Using  the  Lemma  for  x=0,  it  is  clear  that  t  is  unique. 

Then  from  the  analyticity  of  G(0,y)  for  |v|<  1  it  follows  that 

,-10  0  1 
c  Poi'Poi'Poi 

G  (0,0) - r1  ■■■■■  G(0,0)  .  (A15) 

PlO+tPlO 

Substituting  (A15)  in  (A12) ,  G(0,y)  is  determined  up  to  the  constant  G(0,0). 

To  find  G(x,0)  let  f=f(x),  be  the  solution  for  |x|  <  1  of 

e(x,f)  =  0  (A16) 

in  the  circle  |f|  =1.  In  the  lemma  it  was  proved  that  such  a  solution  f  is 


unique. 
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Then  again  from  the  analyticity  of  G(x,y)  for  |xj  <  1  it  follows  that 


Gfx  01  -  c(x,f )G(Q  ,f~)  +d(x , f)G(0 ,0) 

b(x7f) 

Substituting  (A12)  in  (A17) ,  the  function  G(x,0) 
G[0,0).  To  find  G(0,0)  let  x+1  in  (A9)  and 
tion  condition  G(l,l)  =  1  to  obtain"  : 


(A17) 

is  determined  up  to  the  constant 
y-*-l  in  (A10)  and  use  the  normaliza- 


e1(l,l)-b1(l.l)G(l,0)+c1U,l)G(0,l)+d1(l,l)G(0,0)  (A18) 

and 

e2  ^ 1 . 1) =b2 (1 , 1)G(1 ,0) +c2 (1 , 1)G(0, 1) +d2 (1 , 1) G(0,0)  (A19) 


and  from  (A12)  we  have 


G(0, 1)  =  F (0 , 1) 


p^GCO.OMp^+p^G^O.Q) 

1  -  FCO.l)^ 


(A20) 


Solving  (A15) ,  (A18) , (A19)  and  (A20) ,  we  obtain  G(0,0).  The  ergodicity  condition 
is  that  G(0,0)  >  0. 
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Two-node  networks,  (a)  System  1:  non- symmetric  ALOHA  network. 

(b)  System  2:  tandem  network;  no  interference  at  the  station. 

(c)  System  3:  tandem  network;  interference  at  the  station. 
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Fig.  4  -  System  1:  T  versus  p. 


System  1:  p*  versus 
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^min  versus  Y  for  non -symmetric  and  symmetric  access 
schemes. 
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10  -  System  3:  versus  p. 


TWO  NODES 
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12  -  Comparison  of  analysis  and  approximation  of  two  node 
symmetric  ALOHA  network. 


Fig.  14  -  Comparison  of  approximation  and  simulation  results  of 
four-node  symmetric  ALOHA  network. 


